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ABSTRACT 


The results of measurements of thermoacoustically generated 
temperature gradients in short, thin plates located in a resonant tube 
are reported. The temperature gradient results from a heat flow 
generated by the acoustic field. Measurements, at both axial and 
transverse positions in the tube, are made for pressure amplitudes 
ranging from approximately 145 dB to 162 dB re 20 uPa (rms), in 
argon and helium; for single plates and stacks of up to five plates 
having separations from approximately 4 to 40 thermal penetration 
depths. The results are compared to predictions based upon a theory 
by Wheatley and others (J. Wheatley, et al., Journal of the Acoustical 
Society of America, Vol. 74, pp. 153-170, 1983). For pressure 
amplitudes below 150 dB, the measurements agree well with theory. 
At higher pressure amplitudes, the agreement diminishes. It is 
concluded that non-thermal effects are responsible for the 


discrepancies. 


TABLE OF CONTENTS 


DNA tsk) DEP LOC ER VIN eterna cre eae aesesesscessassccsesscecnenecossdenscceacccsnsconsesois 8 

TEER OGG cac gg ten ttre ee enema TOM ETTETTESE NCC Ts cons ososcshesecascocssosccseccsacnsnecceassacssenoreess 10 

ENO aL AIRE Sex@) PO) AG) CIE GOIN Mereenreanet tne neers eve ccc.csiescdeescssscsoessacedeccocedausseesdabccasecs ass 10 

fem rE dec y 18) 2 CO) La) oy TDL OT Oar ee ctineidan J. « ca'vnadialoass «= 2aebtWeaiausivaucasa. se 10 

Pe es clo COMMUN N= (Oa Meme eee 2 Macs duno; cavedsasesceaenoscusdaaccosesssvevonsessess oO 

Ze SUEOUMCINe GaAs COMGICEFAUONS ......ssicessecsevsscacsssscosseessevane 18 

Cee ninenatUnim DIPPMRENCH PREDICTION ...6.0c..00..0cc000se0see0s 19 

AOFM Owe eA PRIVEP NI FEN I @ SSDS Wccrrecsoeoscssvectsonesosssscsesseasescnsessocccccceees 25 

Me EMR 05.605 c2cto sss 4scnseveshneses+cesseaseessusssseiedaseasccaccascccssnsccsesdevecessiieccarssece 27 

fmm ©) Or M Oe Neem VILA EPPING ec. cceacecetece titties: cuss ace cnteseccsssavsectectececsocses 27 

fe ey IDE @ re VIE TING: fics canes. cncscesvactassoansssassecsessacencesenscooessavavscoen 38 

ESRD NE, Ws aces epee es PNP Dlcvassssccossaveeadtcesdoavensescsuueoiies 41 

oe Oe ee Diao PING seeie:....02eete.... 222s... oes sceae es 41 

ese str aly herd aca enV WN MING) 222 ecseeecs ccc csaccesasscessesoscasecetesassesssaseocaseeasacsis 93 

Be a eras OMAN yD lesa rence ee rece tae aaah seu ccousasseuaccsnacesadieoscreyssesesheticeseccssonseen 61 

VI. CONCLUSIONS AND RECOMMENDATIONS ...............cccecececeeeeenees 70 
APPENDIX: PARTIAL LISTING OF THE PHYSICAL PROPERTIES 

Oy Rte) NBN SPE TEU NCI OVI! oo 03 wiiseacecancs svsnsaartsdeusscaasevaveecensouvhsesss Wl 


ACKNOWLEDGMENT 


I would like to thank my thesis advisors, Dr. Anthony Atchley and 
Dr. Thomas Hofler, for their assistance, words of wisdom, and infinite 
patience. Their gentle prodding made this subject interesting and 
much easier to understand. Without their capable hands, my ten 
thumbs and I would still be in the basic construction phase of the 
experiment. 

Next, I thank Mr. George Jaksha of the physics department 
machine shop. His unparalleled ability and insight into sometimes 
poorly worded construction requests resulted in precisely made appa- 
ratus which enabled me to conduct my experiment. 

Next, I thank my fellow students who endured many hours of high 
dBs while trying to study. 

Finally, I thank my wife, Candyce. Her typing and editing skills 
allowed for timely completion of a readable thesis. More importantly, 
however, was her understanding of the many hours required for this 


project and her support through it all. 


I. INTRODUCTION 


Thermoacoustic heat transport has been studied in detail by the 
late John Wheatley and associates [Refs. 1 and 2]. One instrument 
developed in these studies is the thermoacoustic couple, or TAC. 
Wheatley found that, when a TAC is placed in a resonant tube, driven 
so that the acoustic wavelength is much greater than the length of the 
TAC, a temperature gradient is induced across it. The gradient is a 
function of the acoustic pressure amplitude, the type of gas in the 
tube, the thermodynamic properties of the TAC plate material, the 
length of the plate, and the plate configuration (i.e., single or multiple 
plates). 

It is interesting that, while Wheatley states that the thermoacous- 
tic phenomena can be observed in a single plate [Ref. 1:p. 155], his 
single plate results indicate a large deviation from theory [Ref. 1:p. 
158]. This discrepancy is noted, but there is no report of further 
investigation into its cause. 

The purpose of this thesis is to study the thermoacoustic effect, 
on a single plate and on multiple plates, to better understand both the 
basic theory of thermoacoustic heat transport and the reason for single 
plate discrepancies. 

The results of two experiments are reported. In the first experi- 
ment, changes in the temperature difference across a TAC are studied 


as a function of the longitudinal position of the TAC in the resonant 


tube, the number of plates in a stack, the spacing between plates, the 
acoustic pressure amplitude, and the type of gas in the tube. In the 
second experiment, the developed temperature difference is mea- 
sured as the TAC position is varied along a diameter of the tube. The 
results of these experiments are compared to the theory as developed 
by Wheatley and associates. Possible causes for the observed deviations 


from theory are offered. 


Il. THEORY 


A. INTRODUCTION 

A synopsis of the theory of thermoacoustic heat transport, based 
upon John Wheatley’s development, follows. In addition to the basic 
theory, modifications to account for the thermal conductivity of the 
gas will be considered. The possible effects that non-thermal pro- 
cesses may have on the heat transport will also be discussed. The 
reader is directed to References 1 and 2 for a more detailed develop- 


ment of the basic theory. 


B. THERMOACOUSTIC EFFECT 
l. Basic Theory 

In order to understand the thermoacoustic effect, it is neces- 
sary to describe what happens to a gas parcel in an acoustic standing 
wave developed within a resonant tube. In this study, we are consid- 
ering the heat pumping action of such an acoustic wave. 

The four stages of the heat pumping cycle are depicted in 
Figure 1 [Ref. 2:p. 13]. Initially, let us assume that the plate, or TAC, 
has a small temperature gradient, VT, across it. Also, let us assume 
that the idealized plate does not conduct heat in the longitudinal x- 
direction. At step 1, a gas parcel undergoes an adiabatic compression 
from pressure p to p+p,. At the same time, the parcel is displaced 


toward the pressure antinode by a distance x,. This distance may be 
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expressed as xj =p;}/Pmc@ where p} is the acoustic pressure or adia- 
batic pressure change, Pm is the density of the gas, c is the speed of 
sound in the gas, and @ is the angular frequency. Also, at this time, 
the parcel’s temperature increases from T to T+T,. Ty, is the adiabatic 
temperature change and will be discussed in more detail in the fol- 
lowing pages. At step 2, there is a heat flow, Q, in the transverse 
direction, from the hotter gas volume to the cooler plate. This heat 
flow also results in a temperature reduction oT in the gas volume. On 
the return leg, step 3, the gas parcel expands adiabatically and moves 
back a distance x; but is now at a lower temperature T-dT. Therefore, 
at step 4, there is another heat exchange, Q', but this time from the 
hotter plate to the cooler gas. This heat exchange brings the gas back 
to the original temperature T while, at the same time, it reduces the 
temperature of this end of the plate. This heat pumping cycle repeats 
with the frequency of the acoustic wave. 

In describing the heat flow in short stacks associated with 
this process, Wheatley, et al., estimate the hydrodynamically trans- 


ported heat Q" as 


Q" ~ II dx U1 Cp Pm OT (1) 


[Ref 2:p. 14], where 


II = overall perimeter of the plate (measured transverse to the 
flow); 

dy, = thermal penetration depth; 

u; = acoustic velocity magnitude; 

Cp = specific heat at constant pressure of the gas; and 
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oT = change in temperature in the gas due to the heat flow process. 


From Figure 1, step 2, it can be seen that 


omr= 1, —- x) VT =T} (1 Tx (2) 


T,, as earlier stated, is the adiabatic temperature change, and may be 


expressed by 


Tm B 





i 3 
1 = Bmep P! (3) 
where 
Tm = mean absolute temperature; and 
8 = isobaric expansion coefficient. 


For an ideal gas, the product Tmf equals unity. 
The value for VT which makes 6T = O, in equation 2, is 
known as the critical gradient so that 


ae 
VTcrit = ay (4) 


[Ref. 2:p. 14]. This temperature gradient is the largest that can be 
acoustically induced on this plate. 


Using Wheatley’s definition [ = VT/VTcrit [Ref. 2:p. 14] and 
rearranging the terms in equation (1), we can express the hydrody- 


namic heat flow as 


Q" ~ Idx (TmB) pi ui (1 - 1). (3) 
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For the assumed ideal plate, we will remain in the heat 
pumping regime while the heat flow is towards the nearest pressure 
antinode. This will occur as long as VT < VTerit or I < 1. If VT = VTerit, 
then there is no heat flow. 

Now, let us consider that the TAC is not ideal and, thus, does 
conduct heat. There now exists a return heat flow Q,, due to diffusive 
thermal conduction through the plate material. This heat flow may be 
expressed as Q; = kK AVT, where x is the thermal conductivity in the 
direction of the plate length, A is the area of the plate transverse to 
the flow, and VT is the temperature gradient across the plate. The 


net heat flow of this system, Qt, may be expressed as 


Qt = QO" - Gr, (7) 


Or 
Qt = [H 5x (Tm8) piui (1 -1)] -[k A VT]. (8) 


Considering this heat flow, the heat pumping regime is main- 
tained as long as Q" is directed toward the nearest pressure antinode. 
At Q" = Qr, we have reached steady state and Qt = 0. We have, then, 
the maximum temperature gradient, VT, of this system. For a given 
transverse area A, a non-zero x will limit the maximum value of the 
temperature gradient to that achieved when Q" = Qy. 

Now, consider the term II dx (Tm) as constant for any posi- 
tion of the TAC in the standing wave. This results in Q" being depen- 
dent on the p,u, (1 -I) product. [For small acoustic amplitudes, VT 
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is small, = 0, and thus (1 -I) > 1 and is inconsequential. At large 
amplitudes, VT approaches VTe;it and (1 -T) > 0. This will affect Q". 
These two cases will be discussed in more detail later in this paper.] 
By studying equation 8, we can see that, for a given x and plate geom- 
etry, the temperature gradient that exists across the TAC, in the heat 
pumping case, is dependent on the p,u, product. 

The acoustic amplitude p, and particle velocity u, vary with 


position in the standing wave. They may be expressed as 


Pi = Po cos (kx) (9) 
and 
lige g@escim (KX), (10) 


where k is the wave number, x is the displacement from the closed 
end of the tube, and Po and Uo are the magnitudes of the pressure and 
velocity at their respective antinodes. Thus the product p,u, may be 


expressed as 
Po U 
p1uj = | Po Uo) sin (2x) | (11) 


The temperature gradient, therefore, can be expected to vary 
sinusoidally. Figure 2, although not to scale, shows the relationship 
among pj, uj, and VT plus the variation of VTcrit. As is evident, the 
temperature gradient swings from positive to negative values with 
zeroes at either the pressure or velocity nodes. This brings us to the 


point that the heat flow is towards the nearest pressure antinode. 
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This is expected because, as a parcel of gas is moved, in the positive 
x-direction, toward a pressure antinode (Figure 1), the parcel is com- 
pressionally warmed and thus will transfer heat to the plate. If the 
parcel is moved toward a pressure node, it experiences cooling by 
expansion and takes heat from the plate. Furthermore, from equation 
(8) and Figure 2, we deduce that as VTert — O the temperature gradi- 
ent itself must approach zero. 

In our investigation, the plate and gas will initially be at the 
Same temperature. At the onset of the acoustic pressure IT = O and, 
therefore, as seen by equation (5), the heat flow is initially a maximum. 
This will result in a temperature gradient across the plate which will 
increase toward VTcrit. However, as discussed in equations (7) and 
(8), VT will increase only until steady state is achieved, i.e., Qr = Q". 

In Ref. 1 (pp. 155-157), Wheatley, et al., discuss the afore- 
mentioned heat flow in a slightly different manner. There, the devel- 
opment is continued in order to derive an expression with which to 
calculate the temperature difference developed across a TAC. That 


expression is: 
(.25P2& (l+o7) 
(pme (xdg/L) (1 +0) iRaoy nce 9)) sin (2wx/c) 


(1 + Ce dx e -o%)) 
(4(kd2/LipmLTmo (y- 1) (1 - 02) [i= cos(2wx/c)] 


AT = rz) 






ref l:p: 157, Eqn. 17), where 


ava} 
Po 


temperature difference across the plate 
peak acoustic pressure of the driver 
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dx = thermal penetration depth 
o = Prandtl number of the gas 
y = the ratio Cp/Cy of the gas 
L = length of the plate 

Tm = mean ambient temperature 
Om = mean density of the gas 


c = speed of sound in the gas 

® = angular frequency 

Xx = mean position of TAC center measured from the closed end 
of the tube 

« = thermal conductivity of the plate 


dg = thickness of the plate. 
It is this equation and predictions based on it that will be 


examined in this thesis. 


2. Surrounding Gas Considerations 


In the preceding equation, the product xdg2 ignores the 
thermal return path via the gas surrounding the TAC. This return path 
is generally the gas between two plates. Therefore, the thickness of 
the gas path is equal to the plate separation. With a stack of plates, 
this path is well defined and can be expressed as either the absolute 
distance between the plates or as a multiple of the thermal penetra- 
tion depth. Thus, we can express the effective kd2 for a stack of plates 


aS: 


Kdg = (kg * separation) + (Kp * dg), (13a) 


or 


Kd9 = (Kg * dx a5 N) + (Kp = do), (13b) 
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where N is the number of thermal penetration depths equivalent to 
the required separation and kg and kp are the thermal conductivities of 
the gas and plate, respectively. 

For a single plate, the return path is not as well defined. In 
order to estimate the effect of the thermal conduction through the 
gas, we consider two extreme cases. In the first case, the return path 
is assumed to be 5 thermal penetration depths above and below the 
plate, for a total of 10 penetration depths. In the second case, it is 
assumed that the nearest boundary is approximately one tube radius 
away and thus the effective thickness of the return path is one tube 
diameter. In other words, the majority of the space in the tube pro- 
vides the return path. Using this assumption, the effective «dg 


becomes 


Kd2 = (kg * Diameter) + (kp * dg). (14) 


To avoid confusion between the two methods of including the gas 
return path, equations (13a) and (13b) will be referred to as Theory, 
while equation (14) will be referred to as Modified or Mod Theory. It 
is predicted that Theory underestimates the effect of the gas thermal 
conduction, while Mod Theory presents the absolute upper limit on 


the gas conduction. 


C. TEMPERATURE DIFFERENCE PREDICTION 
Now consider a single steel plate TAC, of length 2.26 cm, width 
2.16 cm, and thickness 0.127 mm, placed in an Argon-filled tube at 
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294 K. The thermal conductivity of steel is 17 W/m-K [Ref. 3:p. 
D-187]. The tube is driven at 705 Hz at an acoustic pressure level of 
150 dB re 20uPa, or approximately 600 Pa (rms). The gas return path 
thickness is assumed to be 5 thermal penetration depths for Theory 
and 6 cm for Mod Theory. Using these parameters and the properties 
of Argon [Ref. 4; see Appendix 1], the predicted temperature differ- 
ence can be calculated from equation (12). As can be seen in Figure 3, 
Theory predicts a higher temperature difference than Modified 
Theory. This difference is expected since Mod Theory assumes a 
much larger volume of gas with which to conduct heat. 

Consider now a stack of plates, with plate separation of 1 mm, 
placed in the same resonant tube filled with Argon at 294 K. This 
time, the acoustic pressures will be 150, 162, 168, and 174 dB re 
20uPa. The predicted temperature differences are shown in Figure 4. 
As can be seen, the temperature difference developed across the TAC 
increases with increasing acoustic amplitude. At higher amplitudes, 
there is a tilting or skewing of the temperature curve toward the 
pressure antinode. In contrast, Figure 5 shows the effect of using 
Helium in the tube. The physical properties of this gas [Ref. 5:pp. 48- 
49] result in a much larger temperature difference at the same ampli- 
tudes and also lead to a much more noticeable skewing effect. 

Returning to the discussion on heat flow, one can determine the 
origin of this skewing. In equation (8) and Figure 2, the effect of the 
product piu; was seen. Near a pressure antinode, the pressure pj, is 


large and thus gas parcels undergo a large temperature change due to 
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compression. At the same time, the displacement due to u is small. 
This combination of a large temperature change over a small dis- 
placement means, by definition, that there is a large critical tempera- 
ture gradient. This is confirmed by the large value of VTcrit near the 
pressure antinodes in Figure 2. This also explains the large tempera- 
ture difference near the pressure antinodes observed in Figures 4 and 
5. A large VTerit would allow VT to increase via the (1 - I) term in 
equation (8), and thus for a given TAC length the temperature differ- 
ence would increase. Now, as the TAC is moved from the pressure 
antinode, the value of VTcrit decreases rapidly. This causes a rapid 
decrease in Q", because of (1 - IT) — O, and a corresponding rapid 
change in the developed temperature gradient is predicted. As a 
pressure node is approached, VTcrit approaches zero, but at a slow 
rate. Therefore, a much more gradual change in the temperature gra- 
dient is predicted. This means a corresponding gradual change in the 
temperature difference. Thus, we see that at higher pressure 
amplitudes we expect to see the skewing evident in Figures 4 and 5. 
By considering the separation of the plates as a means to define a 
return path through the gas, it is evident that the predicted tempera- 
ture difference will vary with plate separation. Furthermore, while it 
is not readily apparent from equation 12, Wheatley’s group studied the 
rate of change in temperature on a stack of plates as it relates to plate 
separation distance d, compared to the viscous penetration depth 
distance 5,. Their conclusion was that, “... qualitative changes from 


the simple boundary layer theory may begin to show up for d < 28,.” 
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[Ref. l:p. 155] We would, therefore, expect deviation from theory as 


the separation approaches the viscous penetration depth. 


D. NON-THERMAL PROCESSES 

In addition to the aforementioned thermal effects, there are other 
hydrodynamic effects due to the impinging acoustic wave. While these 
additional processes may have secondary thermal effects, they are 
non-thermal in origin and therefore are identified as non-thermal 
processes. In a standing wave tube which is driven at resonance, 
there is a streaming motion at the tube boundaries which is toward 
the velocity nodes as is depicted in Figure 6 [Ref. 6:p. 348]. As dis- 
cussed earlier, the heat flow developed in this same tube is towards 
the same velocity nodes. If the TAC is on the longitudinal axis of the 
resonant tube, the acoustic streaming will provide an increase in the 
return heat flow. For a single plate, this additional return flow may 
cause the temperature gradient to be less. 

For a stack of plates, the effects of streaming will depend on the 
effect that the additional plates have on the streaming pattern in the 
tube. The plate separation may also be important to the tube pattern. 
The one-way streaming that occurs over the individual plates may be 
more prominent than the tube effects. At this point, we can only say 
that acoustic streaming effects are difficult to predict but they are 
likely to occur at the pressure amplitudes typically used in the investi- 


gation of the thermoacoustic effect. 
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Ill, EXPERIMENTS 


A. LONGITUDINAL MAPPING 

In order to study the temperature variation with respect to posi- 
tion in a standing wave, a copper tube 89.7 cm in length and 6.22 cm 
in diameter was uSed as the resonant cavity. The tube was filled with 
Argon and pressurized to just above atmospheric pressure to ensure 
that all leaks would be outward. Figure 7 shows the equipment 
arrangement for this portion of the experimentation. 

Before data acquisition began, the endplate and probe micro- 
phones were calibrated. The microphone at the endplate (Endevco, 
Model #8550M1) was used to accurately monitor the frequency and 
amplitude of sound in the tube. This microphone is readily calibrated 
using a D.C. calibration method (see Figure 8). The response of the 
microphone was assumed to be flat from D.C. to the experimental fre- 
quencies because the fundamental microphone resonance is at 50 
KHz. The other microphone used was an earpiece from a portable 
television set. As this transducer was used only to monitor the posi- 
tion of the pressure extremes in the tube and not to determine abso- 
lute pressures, precise calibration was not required. Therefore, a 
rough comparison calibration using the Endevco as a standard was 
completed. Later, due to poor performance at high acoustic pres- 
sures, the earpiece was replaced by a similarly calibrated, locally con- 


structed, piezoceramic transducer. 
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To determine a suitable driving frequency, the resonances of the 
Argon-filled tube were studied (see Figures 9 and 10). An interesting 
point of this procedure was the discovery that the speed of sound in 
Argon of 319 m/s at 0° C, obtained from a well-used source [Ref. 3:p. 
E-47], is in error. By comparing the tube length to the observed res- 
onant frequencies, the speed of sound was calculated to be near 320 
m/s at 20° C. Further research found three more-recent sources [Ref. 
4; Ref. 7:p. 3133; Ref. 8:p. 264] which give the value of 319 m/s at 20° 
C. This is the value used for the following experiments. 

Using a single sheet glass TAC, trial runs were conducted to verify 
experimental procedures. This trial provided a method of determin- 
ing the time required for stabilization of the TAC temperature after it 
had been moved (see Figure 11 for typical plot). This was required for 
the automation of data acquisition. Using the measured temperature 
gradient results from this trial run (Figure 12), a decision was made 
regarding the portion of the tube from which useful data could be 
obtained. 

A TAC, consisting of one 0.0127 cm thick stainless steel sheet 
laminated to one 0.0127 cm thick G-10 fiberglass sheet, 2.26 cm long 
and 2.16 cm wide, was fitted with a thermopile consisting of four pairs 
of constantan-chromel thermocouples. To ensure correct voltage 
readings due to the temperature difference developed across the TAC, 
the leads from the TAC are put in thermal contact with a copper plate 


to provide a stable temperature reference (See Figure 7). This “heat 
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sink” arrangement also limited random voltages caused by tempera- 
ture variations along the leads, but not at the TAC. 

The temperature difference developed across the TAC, as a func- 
tion of position from the rigid end of the tube, was measured for 
acoustic pressure levels of 145, 150, and 162 dB re 20 up Pa (rms). 
Next, stacks of 3 and 5 plates were constructed by stacking one or two 
steel plates above and below the measurement TAC at a separation 
distance of approximately 10 thermal penetration depths. Again, the 
temperature difference as a function of position was measured at 
pressure levels 150 and 162 dB. The measurements were then 
repeated for a Separation distance of 36 thermal penetration depths 
and then for a stack of 3 plates separated by 4 thermal penetration 
Sie ptiis. 

To ensure that the data were reproducible, the initial data mea- 
surements were taken twice. Typical agreement in measured values is 
depicted in Figure 13. As more data runs were completed, the repro- 
ducibility of the data was verified. Therefore, data from later runs 
were verified by retaking only random data points. 

The final measurements in this phase of the study were made in a 
Helium-filled tube. The new resonances of the tube were found (see 
Figure 14) and the driving frequency 2120 Hz was chosen to maintain 
the wavelength to TAC length ratio used in the Argon-filled tube. In 
this portion of the experiment, a stack of 3 plates separated by 6 
thermal penetration depths was used for pressure levels of 150 and 


162 dB. 
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B. DIAMETRICAL MAPPING 

Before probing the tube diameter, the results of phase 1 were 
perused to determine a position within the tube where a reasonable 
temperature difference was established at all pressure levels studied. 
This position was 41.5 cm from the endplate, where 2kx = 11.5. The 
tube was then cut at this position, the apparatus shown in Figure 15 
inserted, and the tube sections rejoined. Again, the tube was filled 
with Argon and the resonances were found (see Figure 16). The reso- 
nance frequency closest to the driving frequency used in the first 
experiment was chosen as the driving frequency for this experiment. 
For acoustic pressure levels of 150 and 162 dB, the temperature dif- 
ference across the single TAC was measured as a function of position 
along the diameter of the tube. Then the same measurement was 
made for the center plate of a stack of 3 plates separated by 10 ther- 


mal penetration depths. 
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IV. RESULTS 


A. LONGITUDINAL MAPPING 

As previously described, a TAC constructed of a single glass plate 
was used in a trial measurement. As can be seen in Figure 17, at an 
acoustic amplitude near 150 dB re 20uPa (rms), the measured tem- 
perature difference and acoustic amplitude have the predicted rela- 
tionship {i.e., the temperature difference varies approximately 
sinusoidally with zero crossings corresponding to acoustic pressure 
maxima and minima). 

The temperature difference across the single stainless steel TAC 
was measured as a function of tube position and pressure amplitudes. 
The results are shown in Figures 18 through 20. In these figures, the 
solid line represents the experimental values; the dashed and dotted 
lines represent Theory and Mod Theory, respectively. At acoustic 
amplitudes of 145 and 150 dB re 20uPa (rms), the measured temper- 
ature difference falls between the values predicted by both Theory and 
Mod Theory calculation methods (Figure 18 and 19). At the higher 
level of 162 dB, the peak temperature difference is at least 2° C (or 
approximately 50 percent) lower than that predicted by either 
theoretical method (Figure 20). However, there is still good 
agreement near the pressure antinode (2kx = 12.7). Additionally, 
skewing of the data toward the pressure antinodes is evident, yet not 


predicted. 
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When stacks of 3 and 5 plates, with a plate separation of 36 ther- 
mal penetration depths, are used at an acoustic pressure of 150 dB, 
there again is close agreement between measured and predicted data 
(Figure 21). The agreement between the observed temperature dif- 
ferences for a 3 or 5 plate stack is within experimental fluctuation 
observed on repeat data runs, as discussed earlier (Figure 13). At the 
higher sound level of 162 dB, a large discrepancy, 2° C (about 50 per- 
cent), between measured and predicted peak temperature differences 
is observed (Figure 22), although the agreement near the pressure 
antinode is still good. Again, the 3 and 5 plate stacks produce compa- 
rable results. 

When measurements are made using stacks with plate separations 
of 10 thermal penetration depths, there is again agreement between 
measured and predicted data at the low acoustic amplitude (Figure 23) 
and disagreement at the high amplitude (Figure 24). In this case, the 
magnitude of the disagreement is smaller (only 1° C or about 30 per- 
cent) than the predicted peak. This disagreement is not as much as 
that for a single plate or for a stack with 36 penetration depth spacing. 
Additionally, there again appears to be a strong similarity between the 
results for the 3 plate and 5 plate stacks. 

With a stack consisting of 3 plates, separated by 4 thermal pene- 
tration depths, again there is good agreement between observed and 
theoretically calculated temperature differences at a low level of 
acoustic pressure (Figure 25). Figure 26 shows the now-familiar dis- 


crepancy at the higher driving amplitude, but it has increased back to 
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approximately 2° C (about 50 percent). The noticeable skewing is 
again present, with good agreement near the pressure antinode. 
Finally, this same TAC was used in a Helium-filled tube with a 
plate separation of 6 thermal penetration depths. While the measured 
data are approximately 1° C (50 percent) lower than predicted at the 
150 dB level, the shape of the resulting curve is as predicted (Figure 
27). However, at the 162 dB level, the discrepancy increases to 
approximately 16° C (about 70 percent) lower than the maximum pre- 
dicted. Also, as with the Argon data, the agreement is best near the 
pressure antinode. Furthermore, the observed skewing is much more 
pronounced than predicted (Figure 28). In this figure, the encircled 
jump in the data is the result of a blown fuse in the driver circuit 


which had to be replaced before data acquisition could continue. 


B. DIAMETRICAL MAPPING 

In this experiment, the diameter of the tube was probed with two 
TAC structures. First, using a single steel TAC, the temperature dif- 
ference variations across the diameter were measured for acoustic 
amplitudes of 150 and 162 dB (Figures 29 and 30). At the lower level, 
the temperature difference magnitude increases, albeit less than 
0.01° C (3 percent), as the TAC was moved across the diameter. At the 
162 dB level, the difference magnitude increases about 0.15° C (10 
percent) across the diameter. The probe construction prohibited 
measurement within 0.9 cm of the tube wall. The 3.11 cm distance 


corresponds to the tube center. 
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Next, the measurements were repeated using a 3 plate stack with 
a plate separation of 10 thermal penetration depths (Figures 31 and 
32). At 150 dB, there is a slight decrease in magnitude, 0.05° C (about 
10 percent), across the diameter. There was an interesting result 
observed at 162 dB: The temperature difference went through a 
maximum near the tube center, with a temperature swing of 0.2°C 
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V. DISCUSSION 


The modifications made to Wheatley’s theory allow for considera- 
tion of the return path through the gas, which Wheatley considered 
insignificant for a steel plate TAC. For a single plate TAC at low ampli- 
tude, Theory and Mod Theory are both in good agreement with the 
measured data (Figures 18 and 19). It was stated in Chapter II that 
the Theory and Mod Theory represented lower and upper limits to 
the temperature difference across the TAC. This prediction was con- 
firmed at low pressure amplitudes where the experimental data falls 
between the two predictions (Figures 18 and 19). At the higher 
acoustic amplitude, neither theory can account for the lower tempera- 
ture gradient observed across the TAC (Figure 20). 

It was believed that by studying stacks of plates the cause of this 
discrepancy could be found. For stacks of 3 and 5 plates at various 
separations, the measured data at low acoustic amplitudes compares 
favorably with theory (Figures 21, 23, and 25). This agreement is 
observed no matter the spacing or the number of plates. Even with a 
different gas, there is reasonable agreement (Figure 27). However, at 
the higher acoustic amplitude, the magnitude of the temperature dif- 
ference is noticeably smaller than predicted, regardless of the separa- 
tion, the number of plates, or the gas in the tube (Figures 22, 24, 26, 
and 28). Additionally, the skewing phenomenon is observed. It is, 


therefore, evident that the inclusion of return path through the gas in 
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the theory is not sufficient to account for the discrepancies observed. 
Evidently, at higher amplitudes some other processes are causing the 
discrepancies. 

It is interesting to note the effect that variations in plate separa- 
tion have on the experimental results. As the plate separations are 
reduced from 3 cm, i.e., a single plate, to a separation of 10 thermal 
penetration depths, the agreement with theory improves. But at a 
smaller separation distance (i.e., 4 penetration depths), the agreement 
begins to diverge again (Figure 33). As the separation distance 
decreases from 3 cm to 4 thermal penetration depths, the observed 
temperature difference amplitude peak first increases then decreases. 
Concurrently, the slope of the temperature difference curve first 
steepens, then becomes more shallow. This indicates that there 
exists an optimum plate spacing distance. 

The observed skewing is predicted by the theory (Figures 4 and 
5) but at much higher amplitudes than those used in these experi- 
ments. In addition, theory predicts that the skewing is more pro- 
nounced in Helium than in Argon. Comparison of the results for these 
two gases at 162 dB pressure amplitude reveals that the skewing is 
equally severe. This fact, coupled with the fact that the onset of 
skewing occurs at much lower acoustic amplitudes than predicted, 
indicates that the reasons for the observed skewing are different from 
those presented earlier in this thesis. 

It should be pointed out that the agreement between theory and 


experiment is quite good at all pressure amplitudes in regions near 
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the pressure antinode. In this region, the velocity amplitude is small 
and I is small. The discrepancy is worst in regions of large velocity 
amplitude or large TI (= 1). This observation points to possible regions 
of further investigation in order to explain the discrepancy. 

Theory assumes that the TAC is short compared to the wave- 
length, so that the temperature gradient is constant along the plate 
and that acoustic dissipation due to the plate is small. In our mea- 
surements, the TAC was approximately 5 percent of the wavelength. 
Judging by the close agreement between the results and theory, the 
length of the TAC is sufficiently short for low amplitudes. At higher 
amplitudes, the short stack assumptions may be invalid, resulting in 
the observed discrepancies. 

Still another possible reason for the discrepancies could be that 
the harmonics of the driving frequency are influencing the predicted 
heat flow. Using equation (12) and the acoustic amplitude of the 
second harmonic, measured to be some 20 dB lower than that of the 
resonance frequency, the temperature difference contribution was 
calculated. As can be seen in Figure 34, the maximum temperature 
difference due to the harmonic is only 0.04° C. Therefore, the second 
harmonic is not the cause of the discrepancies. 

If we consider the theory as correct, and if it includes all the 
major thermoacoustic considerations, then some processes of non- 
thermal origin are responsible for the observed discrepancies from 


theory. One such non-thermal process is acoustic streaming. 
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Because acoustic streaming occurs over any surface affected by an 
acoustic wave, we must consider two sources which could influence 
the heat transport. The first source is the streaming associated with 
the tube walls. This flow pattern is as discussed earlier (Figure 6). 
The second source is the streaming generated by the TAC plates 
themselves. It is assumed that the flow generated by the streaming 
would transport heat convectively, similar to what would occur if a 
stream of gas were blown across a plate having a temperature gradient 
across it. If the streaming effects of the tube are dominant, then the 
TAC temperature difference should vary along a diameter of the tube. 
This variation should be manifested as a difference minimum at the 
tube center and at tube walls, with a difference maximum midway 
between them. The measured diameter variations in the temperature 
difference do not demonstrate such behavior. 

The measured data for both a single plate TAC and a stack, at both 
150 and 162 dB acoustic pressure levels, is essentially constant, 
within experimental error, across the diameter (Figures 29-32). The 
magnitude of experimental error was determined by comparing the 
magnitude of the variations to those in successive data runs (Figure 
13), and to the stability of the measured temperature differences 
(Figure 11). It is, therefore, concluded that acoustic streaming gener- 
ated by the tube is not the cause of the discrepancies. Streaming 
effects on the plates occur within a few viscous penetration depths. 
The viscous penetration depth in these measurements is on the order 


of 70 um. The thermocouple wires are 3 mils or 76.2 um in diameter, 
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so the thermocouple joints are larger than the viscous penetration 
depth. Therefore, the effects of plate streaming could not be 
determined. 

The changes in the observed data, when compared to the previ- 
ously measured reSults, indicate an apparent motion of the TAC toward 
the driver. Even when the TAC is purposely positioned 1 cm closer to 
the driver, there still exists an apparent motion (Figures 35 and 36). 
The observed difference in the slopes in Figure 36 is due to the TAC 
starting at positions on opposite sides of the peak of the temperature 
gradient, which exists at the probe position. To determine the cause 
of this apparent motion, the precision of the probe and tube construc- 
tion were studied, as were the variations in driving frequency, driving 
amplitude, and ambient temperature. These changes are considered 
to be of little consequence and the exact reasons for the apparent 


motion are undetermined. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


It is concluded that the theory, for both a single plate and stack 
configurations, is valid at low acoustic amplitudes but not at higher 
amplitudes. Furthermore, assuming that all thermal processes have 
been included in the present theory, it is concluded that non-thermal 
processes, other than acoustic streaming in the resonant tube, are 
responsible for deviations from the theory. Finally, because of the dis- 
crepancies observed, especially at high acoustic amplitudes, it is con- 
cluded that a TAC would be a poor, or at least a very limited, acoustic 
power probe, as suggested by Wheatley [Ref. 2:p. 12], until a better 
understanding of high amplitude effects is achieved. 

To better understand the processes involved and to determine 
the cause of the discrepancies, it is recommended that the experi- 
ments be repeated using: 

a Helium to maximize the thermoacoustic effects; 


b. a shorter TAC to give better agreement with the short plate 
theory; 


c. a phase-lock loop to stabilize the driving frequency; and 


d. various plate separations in a stack to determine the optimum 
plate separation as a function of thermal penetration depth. 
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APPENDIX 


PARTIAL LISTING OF THE PHYSICAL PROPERTIES 
OF ARGON AND HELIUM 
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